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The general relativistic gravitomagnetic clock effect involves a coupling between the orbital motion of a test 
particle and the rotation of the central mass and results in a difference in the proper periods of two counter- 
revolving satellites. It is shown that at 0(c~^) this effect has a simple analogue in the electromagnetic 
case. Moreover, in view of a possible measurement of the clock effect in the gravitational field of the Earth, 
we investigate the influence of some classical perturbing forces of the terrestrial space environment on the 
orbital motion of test bodies along opposite trajectories. 


1 Introduction 

According to Einstein’s relativistic theory of gravitation, the rotation of a mass induces a particular tem¬ 
poral structure in its surroundings that is clearly revealed by the gravitomagnetic clock effect [1-4]. The 
purpose of this paper is to discuss this basic effect and emphasize its close correspondence with an analo¬ 
gous electromagnetic phenomenon, thereby extending our previous work on an alternative derivation of the 
gravitomagnetic clock effect a. Moreover, we discuss the prospects for the observation of this interesting 
effect. 

Gravitomagnetism is expected to be ubiquitous, yet its detection is a very difficult venture. Concerning 
the terrestrial space environment, the GP-B mission a, aimed at the detection of the gravitomagnetic 
precessions of the spins Gl of four gyros carried aboard the spacecraft, has been successfully launched in 
April 2004 and the results are expected to be announced sometime in 2007 111; the anticipated accuracy 
is 1 % or better. A test of another gravitomagnetic effect involving the Lense-Thirring precessions [9-10] 
of the orbits of the LAGEOS and LAGEOS II Earth satellites has been carried out recently [11-12]. The 
claimed accuracy is 5% at one sigma; for a critical discussion of such a test see JTSl . The values of the 
Lense-Thirring precessions of the inner planets of the Solar System have recently been found to be in 
general agreement ITtI with the latest determinations of the planetary perihelia advances ifTSl . although 
the errors are still large. However, they might be measured with a better accuracy in the near future when 
more data will be gathered and processed. 

Among the various gravitomagnetic manifestations which might become observable in the solar system 
within the next few years, the gravitomagnetic clock effect, which involves, among other things, the dif¬ 
ference of the sidereal periods of two free counter-orbiting particles about a rotating massive body, is of 
special interest since it contradicts the common notion of the “dragging of inertial frames”. In the follow¬ 
ing we will show that this effect has a simple correspondence in the electromagnetic case by considering 
the motion of a test charge in a dipolar magnetic field, thereby making the nature of the clock effect more 
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Fig. 1 A test charge (m, q) moves about a central slowly rotating body (M, Q). The figure depicts the 
initial configuration when the rotation of the source is “turned on”. The unperturbed orbit of the test mass is 
assumed to be circular; moreover, ipo is the initial angle in the orbital plane, while (j>o is the corresponding 
angle measured along the equatorial plane to which the orbital plane is tilted by the inclination angle i. 


transparent (see also Q). Finally, the paper ends with a brief discussion of the required accuracy in the 
determination of the orbital elements in order to test the clock effect in a possible experiment with two 
counter-revolving satellites about the Earth. 

2 The electromagnetic scenario 

We consider a spherical body of homogeneously distributed mass M with a charge Q, which rotates with 
a small angular velocity uJq, thereby producing a weak magnetic dipole moment The axis of rotation 
and the direction of shall coincide with the positive z-axis of the background inertial coordinate system. 
We study the slow motion of a test charge of mass m M and charge q {qQ < 0) about the central body 
such that in the absence of the perturbation due to /x^, the motion takes place along a circular inclined 
orbit, where the nodal line coincides with the a:-axis, as shown in Fig. 1. Once the magnetic perturbation 
is ’’turned on”, the orbit will no longer be planar; however, the orbital radius can remain constant to first 
order in the perturbation [5]. Throughout this paper it will be assumed that the test particle moves along 
a spherical orbit of constant radius and that the magnetic field is weak so that second-order effects may 
be neglected. Moreover, we will ignore all radiative and special relativistic effects on the motion of the 
electric charges. 

In the following we show that the motion of the test mass is modified due to the weak dipolar magnetic 
field in two ways: (a) the orbital plane precesses about the spin axis of the central body, and (b) the along- 
track velocity of the test particle is changed. 


2.1 The equations of motion 

The rigid rotation of the central charge produces a current density j which involves a magnetic dipole 
moment 



(1) 
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and hence via the magnetic vector potential 


A = 


J.3 ’ 


the magnetic dipole field 


B = V xA = 


3r(r ■ /Xq) - 


( 2 ) 


( 3 ) 


The motion of the test mass in the electric and magnetic field produced by the central charge is governed 
by the Lorentz force, 


F = q 


E -\— r X B 


( 4 ) 


where r is the position vector joining Q with q, an overdot denotes differentiation with respect to time, 
E = {Q/r^)r is the electric field and B is given by Q. For later purposes, we formally introduce electric 
and magnetic charges Qe, (7e and Qb, (Zb, respectively, and write the equations of motion as 


r = 


Qe^e r 
m 


+ —rxB, 
me 


( 5 ) 


with Q = Qe = Qb and (z = Qe = (Zb- It must be emphasized that the magnetic charge (such as Qb or qg) 
introduced here has no connection with the magnetic monopole strength that is assumed to be zero in this 
work. 


2.2 Precession of the orbital plane 
Since the orbiting test charge possesses a dipole moment 




2mc 


(6) 


where I is the orbital angular momentum, a torque M = /x ^ x B will act upon leading to a precession 
of the orbital plane about the z-axis with angular velocity FJ,, according to £ = SI, x £. 

In order to find $7,, we first determine the vector potential A. Since the magnetic field B is associated 
with the uniform rotation of the central charge of angular momentum J = where I is the moment of 
inertia of M, the dipole moment of the central spherical object becomes 


_ Qb j _ Qb^ 

“ 2Mc ~ 2Mc^' 


or, upon inserting into © 


( 7 ) 


A = aujq X r, 


( 8 ) 


with a = QbI / {2Mcr^). Vector multiplication of (|5]l with mr and using B = aV x (otg x r) as well as 
substituting p = mr in the canonical momentum P — p + {qg/QA yield 


L- 


9b 


(f 


r X A + r X A] = 


29b Q( 


r X (r X tUq), 


( 9 ) 


where L — r x P. Finally, since lJq ^nd r do not depend upon time, we arrive at 


L = —X (r X r) = r2, X L + J x (r x A), 

c ^ ^ ’ 2Mmc^r^ ^ " 


(10) 
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where 


n 




Qua 

- UJ 

me 


Q 


QbQb J 
Mm 2c^r^ 


( 11 ) 


The last term in equation (10) is of since A = and is second order in the magnetic 

perturbation; therefore, it is negligible in accordance with our approximation scheme. It follows that the 
canonical angular momentum precesses with frequency 17, that is “slow” in comparison with the “fast” 
orbital motion. 

We note that the product will always be negative, since we are considering attractive forces only, 
therefore 17, and J point in the same direction. It is interesting to note that in our approximation scheme, 
the orbital speed as well as the magnitude of the orbital angular momentum remains constant in accordance 
with equation (5). What we have demonstrated here for canonical angular momentum also holds true, 
when averaged with respect to the “fast” orbital motion, for £ = r x p. That is, it can be shown that on 
the average, the orbital plane of the test charge is “dragged” in the same sense as the rotation of the central 
charge. 


2.3 The change in the orbital velocity 

In the following we will calculate the change in the along-track velocity due to the magnetic dipole held Q. 
The result can be understood heuristically by determining the magnetic hux through the orbital plane. If 
we assume that the magnetic held is switched on and increases from zero to B, then by means of Faraday’s 
law we can calculate the electric held strength E induced in the orbital plane due to the varying magnetic 
held and obtain 


where i is the inclination angle. From (fT^ and QbE = mdv/dt = mrdip/dt, the change in the angular 
velocity A?/; caused by the dipolar magnetic held becomes 


A'l/) 


QbQb j 

T7 - , r, , cost; 

Mm 


(13) 


this result can be obtained directly starting from equation (5). Since Qb^b < 0, we see that the test charge 
moves slower in the prograde direction (i.e. in the direction of the rotation of the source), and faster in the 
retrograde direction as compared to a situation without the magnetic held. This behavior is also intuitively 
conhrmed by the subsequent evolution of the motion under the action of the Lorentz force: for the prograde 
motion the Lorentz force points away from the source and weakens the attractive Coulomb force. We note 
that this feature is independent of the sign of the revolving charge, because the exchange of the signs of the 
central and orbiting charge will also change the direction of the magnetic held. 

The change in the velocity of the revolving charge can also be anticipated by virtue of Lenz’s rule: the 
induced current will have such a direction as to counteract its own source. Indeed, let us replace the rotating 
source by a single current loop and the orbit of the revolving charge by a conducting wire concentric to 
the inner loop. Any change A/ of the inner current will induce a current in the outer wire having the same 
direction as AJ, thus leading to an increase of the magnetic Hux through the inner loop. This additional 
hux in turn induces a current opposite to A/ and thus opposes a change of the inner current. 


2.4 Time of revolution 

According to (fT3l) . the velocity of the test charge depends on the direction of the motion and thus its 
period will be different for pro- and retrograde orbits. If the magnetic held is absent, the Kepler period 
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is Tk = 27r/'0K, with the frequency t/^K = {\Q-Eq-E\/, while in the presence of a magnetic field, 
'0 = '0K + At/) and the time of revolution T* becomes 


^ ^ ^ ^ A _ 1 gsftQ 

-Ip tpK\ 'fpK ) t/’K V 'Ipf,2mcr^ 


(14) 


The difference in the period for two orbits with different inclinations zi and *2 is immediately found from 
(fT4l) and (|7| to be 


rjn* rj~i^ 


1 J 
2^ Mc^ 


(cosii 


cos 12 ) 


(15) 


for electric charges. Note that this time difference, which refers to the closure of the angle ip in the orbital 
plane, depends only on the specific angular momentum of the source, but neither on the charge nor the 
distance r of the test mass, and it is maximum for opposite equatorial orbits (i.e. zi = 0, t 2 = tt or vice 
versa). 


2.5 The magnetic clock effect 

As demonstrated in the preceding paragraphs, the motion of a test charge in a central Coulomb field to¬ 
gether with a weak dipolar magnetic field results in a slow precession of the orbital plane as well as a 
slightly increased or decreased along-track velocity. 

An inertial observer will define the period of the particle as the time elapsed between two successive 
crossings of a fixed plane in space, and consequently this period will differ from (fl4l i due to the nodal pre¬ 
cession of the orbital plane. The time of revolution as measured by the inertial observer can be calculated 
by noting that the particle’s position angle ip in the orbital plane is related to the corresponding angle p 
measured from the lines of the ascending nodes in the equatorial plane by 


tanp = tan ■)/) cost. 


(16) 


Now, let ipo = ip{t = 0) and pQ = p{t = 0) denote the initial position of the particle and T be its period 
with respect to the inertial frame, where T = T* + AT and AT <C T, T*. The angle swept over by the 
particle after completion of one revolution as seen by the inertial observer in the orbital plane becomes 
po + pT — uipr cos i, while in the equatorial plane this corresponds to pQ -\-2 tt — a;,T and therefore 


tan(0o ~ '^qT) = tan(?/)o + pT — uj^T cosi) cosi. 


(17) 


By means of trigonometric relations and by noting that tan pT ~ pAT, tan a;,T ~ a;,T* and p ^ ^ 
pp^{\ — Ap/pP) together with (fTTTi and (fOT l. we finally obtain to first order by virtue of (fThl l 


27r 

T = ^ 
pK . 


1-^ _ 2tan^icos^ po) cosi 

p^ 2mcr-^ 


(18) 


From this equation, the difference in the periods for two orbits with inclination T and 12 is found to be 


Ti - T 2 = 


TT^b fig 


[2 cos^ po (ta,iP i \cos i\ — tan^ cos ip) — cos p + cos iP) 


|Qe9e| c 

which reduces for = gg and Qe = Qb to the magnetic clock effect 

Ti — T 2 = —-TT „ „ [2 cos^ po (tan^ ii cosii — taiP p cosp) — cosii -f cosi 2 ] ; 

2 Mc^ 


(19) 


(20) 
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for opposite or supplementary orbits {ii = i, i 2 = tt — i), this relation simplifies to 


Tl - To = TT 


Mc2 


(l — 2 tan^ i cos^ V'o) cos 


( 21 ) 


It should be emphasized that the dependence of the period (and of the time difference in turn) on tpQ is due 
to the non-linear relation (fTbl l. 

It is interesting to note that the period for a prograde orbit is longer than the Kepler period for low 
inclinations (remember that qg/ig < 0) and shorter for high inclinations, because in the first case the lower 
orbital velocity implies T > T^, while in the latter case the nodal precession, which makes the period to 
decrease, dominates. 


3 The gravitomagnetic scenario 

In the case of slow motion and weak gravitational fields, the spacetime metric can be approximated by 

P/ll/ — P/ll/ “b ^/ii/i (22) 

where = diag(—1, -|-1, -|-1, -|-1) is the Minkowski metric and |/i/ii/| <C 1 is a small correction so that 
only terms linear in will be of significance. 


3.1 The equations of motion 

As in the electromagnetic case, the mass current density jg associated with the rotation of the central mass 
with angular momentum J produces a dipole moment 


1 /■ • ^3 G r 

y / V ^ 7-' 

with the corresponding vector potential 

. G J X r 
An = - 


9 ^ ^3 


which implies the gravitomagnetic dipole field 

G 3r(r • J) — Jr^ 


Bg =V X Ag = 


(23) 


(24) 


(25) 


analogous to Q. 

In the case of weak perturbations (l22l i and upon restricting to static fields, the geodesic equation can be 
written in the form 0| 


F = —m 


E, 


€g + -r X Bg 

c 


(26) 


where Eg = GMr/r^ and Sg — 0{c~^) is the post-Newtonian gravitoelectric field. As has been first 
pointed out in ifTbl . a complete correspondence with electromagnetism can be achieved by introducing 
the notion of gravitoelectric and gravitomagnetic charges, and q^, respectively, and ignoring £g. The 
equation of motion (l26l l thus becomes 


r = 


GMq^ r 
m 


+ ^rxBg, 
me 


(27) 


in perfect agreement with Q. The appearance of the factor 2 in the “magnetic” force term of (|26] | suggests 
= —2m as well as q^ = —m, where the additional minus sign accounts for the fact that gravitational 
charges of opposite polarity, like in the electromagnetic case, should be attractive. From (l27T i we also read 
Qe = GM and consequently Qb = 2GM. 
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3.2 The gravitomagnetic clock effect 

Because the equations of motion of electrodynamics and linear general relativity are formally equivalent, 
in the weak-field and low-velocity limit gravitational phenomena are expected to be similar to those known 
in electrodynamics. Therefore, by means of the substitutions derived in the preceding section, 

gE ^ -TO, Qe GM, 

Qe ^ -2m, Qe ^ 2GM, 


we can directly translate electromagnetic phenomena into the corresponding gravitoelectromagnetic phe¬ 
nomena. It is important to note that the ratio of the gravitomagnetic charge to the gravitoelectric charge 
is always 2, which corresponds to the spin-2 character of linear gravitation. This is in agreement with 
the fact that the ratio of the magnetic charge to the electric charge is always unity in Maxwell’s spin-1 
electrodynamics. 

With /Tq ^ GJ/c from (|7]i we immediately obtain by virtue of (fTTT i 


^LT 


2GJ 


(29) 


which is the well-known Lense-Thirring precession of the node of the orbital plane of a satellite. Further, 
let T± denote the period of revolution in the same (H-) and opposite (-) sense as the rotation of the source, 
then (fTsT i translates into 


r± = Tk ± 2tt ^ ^ (l ~ 2 tan^ i cos^ ijjo) cos i, (30) 

which is equivalent to equation (35) of ||2l if terms of 0{c~‘^) are neglected. Finally, for opposite equatorial 
orbits it follows from (l30l ) that 


r+ - T_ = dTT 


j 


(31) 


yielding the gravitomagnetic clock effect obtained in m. It is an interesting fact that (ImT i. as well as 
its magnetic counterpart (l2Tli . neither depends on the coupling constant nor on the distance r of the test 
masses. The reason for this is seen by inspection of (fT3T i and (fT4li : both the square of the Kepler frequency 
and the induced velocity change depend on G/r^ (Qe 2GM in (foT l) and hence it will not appear in the 
expression T+ — T_ cx Aijj/tpl. This result is a consequence of the r ^ fall-off of the dipolar magnetic 
field. In a homogeneous magnetic field, for instance, Aip will be independent of the distance and the time 
difference will thus become T+ — T- (x r^. Further, the absence of the gravitational constant G in (l3Tl l 
allows the time difference to become relatively large; inserting the appropriate values for the Earth into 
(EB yields T+ - r_ - 10-^ s. 


4 An alternative clock effect 


If we introduce the angle ^ = w -f cos i + M., where tu, fl, At and i are the argument of pericenter, the 
longitude of the ascending node, the mean anomaly and the orbital inclination, respectively, its temporal 
rate of change is given by 


dt 


du) 

dt 


dfl . di 

H—— cos i — il sin i— + 
dt dt 


dM 

dt 


(32) 


Upon inserting the secular contributions of the oblateness of the Earth as well as the corresponding 
Schwarzschild and Lense-Thirring rates we find that the time needed for ^ to cover an angle of 27r is 
given by 


27r 

n 


I + J 2 


3i?2 (1 -b 
4a^(l — 


(1 — 3cos^ i) — 




4nJ cos i 


0 ^( 1 -e^) c^M(l — 


, (33) 
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where a, e and n are the semimajor axis, eccentricity and mean motion, and J and M are the radius, 
angular momentum and mass of the Earth, respectively; J 2 is the lowest zonal harmonic coefficient in 
the multipolar expansion of the terrestrial gravitational field. For satellites along identical but opposite 
orbits, all but the last term in (l3^ cancel when taking the time difference in the periods of revealing the 
gravitomagnetic signature 


+ _ ^ ^ IGttJ cosi 

« « c2M(l-e2)3/2’ 


(34) 


which could in principle be measured by means of two counter-revolving satellites around the Earth. We 
note that (l3^ reduces to equation (14) of ifTTl in the case of nearly equatorial and circular orbits. 


5 Error analysis 

Equation (1?^ holds only for the ideal case that the two orbits are identical and the orbital elements are 
exactly known. In a real experiment, both the orbits will differ from each other due to injection errors and 
the elements are only determined with limited accuracy. These residuals will contribute to the difference 
in equation (l34t and must hence be kept less than the Lense-Thirring difference. Although equatorial 
orbits would be favorable, however, since most satellite laser ranging stations are located in the northern 
hemisphere, inclined orbits are more suitable from an observational point of view. For a nominal orbit 
with a = 12 000 km, e = 0.01 and i = 63.4°, which represents a frozen perigee configuration, the 
upper panel in Fig. 2 displays the maximum allowable uncertainty of Aa, Ae and Ai, as a function of 
the corresponding orbital elements, while the lower panel illustrates the maximum allowable change of 
A{GM)/GM and AJ 2 /J 2 with respect to the deviation from the nominal orbit of a = 12 000 km. As 
can be seen, a relative accuracy of Aa/a < 6 x 10 ~^ 2 ^ Ae/e < 1 x 10“^, Ai/i < 8 x 10“® is required, 
while the uncertainties in the values of GM and J 2 are of minor importance, since they are known to an 
accuracy of A(G'M)/GM ~ 10“® and AJ 2 /J 2 ~ 10~^ ifTSl . With about 2500 revolutions per year, after 
4 years of observation and for a target error of 1%, the semimajor axis must be known with an accuracy 
of about 1 part in 10®, i.e. to a few cm, e to 1% and i to 0.0001%. The stringent limit for i is due to the 
high inclination of the nominal orbit; for near equatorial orbits this limit becomes less severe Qt). Finally, 
we note that care must also be taken in the accurate modelling of all non-gravitational perturbations. In 



0 200 400 600 800 1000 0 200 400 600 800 1000 

Deviation [m] Deviation [mj 

Fig. 2 Maximum tolerable error of various values for a nominal orbit of o = 12 000 km, e = 0.01 and 
i = 63.4°. 
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this case, however, one would greatly benefit from the experience gained in the LAGEOS-LAGEOS II 
Lense-Thirring experiment. 

6 Conclusion 

It has been shown that the gravitomagnetic clock effect for spherical orbits has a simple correspondence 
with the electromagnetic case. The difference in the sidereal periods of two counter-revolving charges 
about a central spinning charge is due both to the different velocities of the test charges and to the precession 
of the orbital plane. The slower prograde velocity, which in the gravitational case is against the common 
notion of the dragging of inertial frames, can be considered as a consequence of the Lenz rule. The 
gravitomagnetic clock effect is obtained by simply substituting the electric and magnetic charges by the 
corresponding gravitoelectric and gravitomagnetic ones. With regard to the gravitational perturbations, a 
key issue for a possible observation of the clock effect with two counter-orbiting satellites is the stringent 
requirement in the accuracy of the inclinations. 
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